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ABSTRACT 



There is increasing evidence that many galaxies host both a nuclear star cluster (NC) and a super-massive black 
hole (SMBH). Their coexistence is particularly prevalent in spheroids with stellar mass 10 8 - 10 10 M . We study the 
possibility that a stellar-mass black hole (BH) hosted by a NC inspirals and merges with the central SMBH. Due to the 
high stellar density in NCs, extreme mass-ratio inspirals (EMRIs) of BHs onto SMBHs in NCs may be important sources 
of gravitational waves (GWs). We consider sensitivity curves for three different space-based GW laser interferometric 
mission concepts: the Laser Interferometer Space Antenna (LISA), the New Gravitational wave Observatory (NGO) and 
the DECi- hertz Interferometer Gravitational wave Observatory (DECIGO). We predict that, under the most optimistic 
assumptions, LISA and DECIGO will detect up to thousands of EMRIs in NCs per year, while NGO will observe up 
to tens of EMRIs per year. We explore how a number of factors may affect the predicted rates. In particular, if we 
assume that the mass of the SMBH scales with the square of the host spheroid mass in galaxies with NCs, rather than 
a linear scaling, then the event rates are more than a factor of 10 lower for both LISA and NGO, while they are almost 
unaffected in the case of DECIGO. 
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1. Introduction 

Nuclear star clusters (NCs) are massive star clusters located 
at the centre of galaxies (see, e.g., Boker 2010 for a recent 
review). They are common and quite ubiquitous, as the 
fraction of galaxies with unambiguous NC detection is ~ 75 
per cent in late-type (Scd-Sm) spirals (Boker et al. 2002), 
~ 50 per cent in earlier-type (Sa-Sc) spirals (Carollo et al. 
1997) and ~ 70 per cent in low and intermediate luminosity 
early-type galaxies (Graham & Guzman 2003; Cote et al. 
2006). The typical size of NCs is similar to that of Galactic 
globular clusters (GCs, half-light radius w 2 — 5 pc, Geha, 
Guhathakurta & van der Marel 2002; Boker et al. 2004; 
Cote et al. 2006), but NCs are 1 — 2 orders of magnitude 
brighter (and more massive) than Galactic GCs (Walcher 
et al. 2005). Unlike most of the Galactic GCs, NCs have 
often a complex star formation (SF) history, with multiple 
episodes of SF (Rossa et al. 2006; Walcher et al. 2006). 
NCs are of course not the Bahcall-Wolf (Bahcall & Wolf 
1976, 1977) cusps predicted to exist around SMBHs: NC 
sizes are much larger (e.g. Merritt & Szell 2006) and they 
are observed in galaxies thought not to host SMBHs (e.g. 
Ferrarese et al. 2006). 

NCs obey scaling relations with host-galaxy properties 
(such as galaxy mass, spheroid luminosity and velocity 
dispersion, e.g., Balcells et al. 2003; Graham & Guzman 



2003; Balcells, Graham & Peletier 2007; Ferrarese et al. 
2006; Wehner & Harris 2006; Rossa et al. 2006; Graham & 
Driver 2007). These scaling relations are similar to those 
observed for super-massive black holes (SMBHs), suggest- 
ing that there is a link between NCs and SMBHs. Various 
studies (Ferrarese et al. 2006; Wehner & Harris 2006) in- 
dicate that SMBHs (NCs) are found predominantly above 
(below) a stellar mass threshold M* ~ 1O 1O M . However, 
the presence of a SMBH and that of a NC do not seem 
to be mutually exclusive. In addition to the Milky Way 
(MW, see, e.g., Schodel et al. 2007, 2009), Filippenko & 
Ho (2003) identified at least one galaxy (NGC 4395) host- 
ing both a SMBH and a NC, and Graham & Driver (2007) 
subsequently reported the existence of two additional such 
galaxies (NGC 3384 and NGC 7457). The sample of galaxies 
hosting both a SMBH and a NC was substantially increased 
by Gonzalez Delgado et al. (2008, 2009), Seth et al. (2008a) 
and Graham & Spitler (2009). In particular, Graham & 
Spitler (2009) identify a dozen NCs in galaxies that host 
SMBHs, whose mass was already determined via dynamical 
measurements. Furthermore, Graham & Spitler (2009) sug- 
gest that most spheroids with stellar masses ranging from 
10 8 M Q to a few 10 10 M o might host both a SMBH and a 
NC. 

NCs have been proposed as important sources of gravi- 
tational waves (GWs, Poincare 1905; Einstein 1916, 1918), 
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originating from mergers between stellar-mass black hole 
(BH) binaries (Miller & Lauburg 2009; see also Freitag 
2003; Hopman, Freitag & Larson 2007; O'Leary, Kocsis & 
Loeb 2009; O'Shaughnessy, Kalogera & Belczynski 2010). 
In fact, because of their relatively high escape velocity, NCs 
should be able to retain most of their stellar-mass BHs 
(Miller & Lauburg 2009), in spite of natal kicks and gravita- 
tional interactions (although the presence of a SMBH inside 
the NC might affect this scenario, see, e.g., Merritt 2006, 
2009). Furthermore, three-body interactions are expected 
to be very frequent in NCs, given their high central den- 
sity, leading to shrinking and coalescence of hard binaries. 
For these reasons, Miller & Lauburg (2009) predict that 
tens of BH— BH mergers per year in NCs will be observable 
with the advanced Laser Interferometer Gravitational- wave 
Observatory (advanced LIGO, Harry et al. 2010). Given the 
increasing evidence of the coexistence between SMBHs and 
NCs in some galaxies, we propose that mergers between 
SMBHs and stellar BHs may also occur in NCs. Such events 
fall within the class of extreme mass-ratio inspirals (EMRIs, 
see Amaro-Seoane et al. 2007, for a review) and may be de- 
tected by space-borne interferometers. At present, it is still 
fairly uncertain what mission could fly in the future. As ref- 
erence, we consider the published sensitivity performance 
of the Laser Interferometer Space Antenna (LISA, Bender 
et al. 1998; Gair, Tang & Volonteri 2010), the European 
New Gravitational wave Observatory (NGO, derived from 
the previous LISA proposal, Amaro-Seoane et al. 2012a, 
2012b) and the DECi-hertz Interferometer Gravitational 
wave Observatory (DECIGO, Kawamura et al. 2006). 

In this paper, we investigate the capture rate of BHs 
onto SMBHs in NCs which yield EMRIffl and we estimate 
their detectability by future gravitational- wave space-based 
observatories. 



2. Method 

In this section, we derive the rate of EMRIs in NCs and 
we infer an expected detection rate of such events by LISA, 
NGO and DECIGO. 



2.1. The detection rate 

The detection rate R of EMRIs in NCs can be expressed as 
(see, e.g., Miller 2002; Mapelli et al. 2010) 



R 



dmsMBH x 



L 



2max(wiSMBH,m co ) 



d 3 AWi dt e dV c 



' dmsMBH d< c dVc dt Q dz ' 



(1) 



where m\ and 7712 are the minimum and maximum 
SMBH mass (ttismbh); A^emri is the number of EMRIs 
of SMBH— BH binaries; t e and t Q are the time in the rest 
frame of the source and of the observer, respectively (thus, 



1 In the following, we will use the definition of EMRIs as 
merger events occurring from an orbit with semi-major axis, 
at the moment of capture by the SMBH (the capture radius is 
defined as 8r g , r g being the event horizon of the SMBH), less 
than 10~ 5 pc (Merritt et al. 2011). 
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where c is the light speed, Ho is the Hubble constant (Ho = 
71 km s" 1 Mpc- 1 ) and £{z) = [Q A + (1 + zfVt M ] 1/2 , with 
tt A = 0.73 and D, M = 0.27 (Larson et al. 2011). 

The quantity z max (msMBH, fn co ) is the maximum red- 
shift at which an event can be detected with a sky-location 
and orientation averaged signal-to-noise ratio (SNR) > 10 
by a single interferometer. Thus, z max (msMBH, fn co ) defines 
the instrumental horizon of a given detector. In observa- 
tions with a network of instruments, the signal-to-noise ra- 
tio scales as the square root of the number of instruments, 
and in this respect the results presented here should be con- 
sidered as conservative. The maximum redshift depends on 
(i) the mass of the SMBH tosmbHj (h) the mass of the 
companion (m co ) that merges with the SMBH, (iii) the or- 
bital eccentricity e(t), (iv) the spin parameter of the SMBH 
(j = 5 c/G m| MBH , where S is the modulus of the angular 
momentum of the SMBH, c is the speed of light and G is 
the gravitational constant), and on various other parame- 
ters, as well as on the sensitivity of the instrument. In this 
paper, we adopt m co = 10 M©, corresponding to the typical 
mass of a stellar-mass BH. We consider different values of 
the eccentricity and of the spin parameter. See Appendix A 
for details on the derivation of z max (msMBH, m co ). 

The EMRI rate per unit mass, time and volume can be 
written as: 



d 3 iV EM Ri 



dmsMBH di dV c 



^EMRI 



dn s 



MBH 



dmsMBH ' 



(3) 



where usmbh is the comoving density of SMBHs sur- 
rounded by NCs and i^emri is the rate of EMRIs (i.e., the 
rate of mergers that end-up as EMRIs) per SMBH. 

Graham & Spitler (2009) indicate that SMBHs coexist 
with NCs in spheroids with mass 10 8 < TO sp h/M© < 10 10 . 
The questions (i) whether there is a scaling between the 
mass of the SMBH and that of the host spheroid, and (ii) 
whether such relation is valid for all SMBH masses have 
been debated for a long time. Marconi & Hunt (2003, here- 
after MH03) find that msMBH ~ 0.002 m sp h, but there 
might be large deviations (up to a factor of 10) from this 
value for low-mass (< 10 6 M Q ) SMBHs (see, e.g., tab. 2 
of Graham & Spitler 2009). Graham (2012a, hereafter 
G12a) and Graham (2012b, hereafter G12b) propose that 
the mass of the SMBH scales (nearly) with the square of 
the host spheroid mass in galaxies with NCs (i.e. galax- 
ies with m sp h < 10 11 M©), while it scales linearly with 
the host spheroid mass in high-mass galaxies (which of- 
ten have partially depleted cores). In particular, G12a pro- 
poses that log lo (m S MBH/M ) = (8.38 ± 0.17) + (1.92 ± 
0.38) log 10 (m sph /7 x 10 10 M Q ) for galaxies with NCs. The 
BH mass range corresponding to 10 8 < m sp h/M© < 10 10 
is 2 x 10 5 - 2 x 10 7 Mq and 8 x 10 2 - 5.7 x 10 6 M©, ac- 
cording to the results by MH03 and by G12a, respectively. 
In the following, we will consider both cases. In particu- 
lar, we define / B HK P h) = 0.002 m sph and /BH(m sp h) = 
10 8 38 M Q (m sph /7 x 1O 1O M ) 192 , when adopting the for- 
malism by MH03 and by G12a, respectively. Therefore, 
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equation ([T]) can be rewritten as 

3 /■'"aph2 



i? = 47T | — 

Mo 
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where m sp hi = 10 s M Q and m sp h2 = 10 10 M© are the 
minimum and the maximum mass of spheroids that host 
SMBHs together with NCs, respectively; /bb is the frac- 
tion of spheroids in the considered mass range (10 8 < 
m sph /M Q < 10 10 ) that host both a SMBH and a NC (we 
assume /bb = 1, according to Graham & Spitler 2009); n sp h 
is the comoving density of spheroids. In equation 31 we used 



the fact that , d " SMBH 

amsMBH 



/l 



bb ■ 



dn. 



dmsMBH dmsph ' 

We can derive n sp h from the comoving density of halos 
Uh, assuming that spheroids with mass m sp h are located in 
halos with total mass mi, = m sp h/p. Thus, equation Q can 
be rewritten as 



3 r m h 



m hl 

Zmax [/bh (p m h ) ,m co ] 



1 



dm h ^emri x 



dz /bb x 



dmh 



(l + z)£(z) \J £{z) 



(5) 



where 4^ is the Press-Schechter function (Press & 
Schechter 1974; Eisenstein & Hu 1998, 1999), m hl = 
msphi/p and mh2 = m sp h2/p- We use p = 1.06 x 10 -2 as 
a fiducial value, but we consider also p — 5.04 x 10~ 2 and 
= 2.81 x 10 -3 as upper and lower values, respectively, de- 
pending on the baryonic fraction in dark matter halos (see 
Appendix B for the derivation of p) . In the equation Q , we 
write /bb inside both integrals (in mass and redshift), al- 
though, in the following, we will assume for simplicity that 
/bb is a constant. 



2.2. The EMRI rate 

The EMRI rate ^emri (i-e. the rate of mergers that end-up 
as EMRIs) is likely the most difficult quantity to estimate 
in equation (J3J) , as it depends on the dynamics and on the 
relativistic effects in the neighborhoods ( ^ 10~ 2 pc) of the 
SMBH. In particular, the collisional dynamics of relativis- 
tic star clusters around SMBHs is poorly understood. The 
presence of the SMBH potential well strongly affects the 
inner regions of the star cluster and the dynamics of binary 
stars (Alexander 1999, 2003, 2005; Alexander & Hopman 
2009; Hopman 2009). To derive j^emri in an accurate way 
is beyond the aims of this paper. Here, we will adopt ap- 
proximate rates derived by Merritt et al. (2011). Capture of 
stars on EMRI orbits can be driven by resonant relaxation 
(Rauch & Tremaine 1996; Hopman & Alexander 2006a), 
by two-body relaxation (Merritt et al. 2011) or even by 
dynamical friction (Antonini & Merritt 2012), that is by 
physical processes where background stars exert torques on 



the stellar orbits. The capture of a star is hampered by the 
so-called Schwarzschild barrier (Merritt et al. 2011), that is 
by an angular-momentum barrier associated with the value 
of the orbital angular momentum at which the background 
torques become ineffective due to the Schwarzschild pre- 
cession of the orbifl Merritt et al. (2011) analyze various 
mechanisms to penetrate the Schwarzschild barrier (based 
on classical non-resonant relaxation), and derive a time 
scale ii oss for stars to 'definitely' cross the barrier (see their 
equation 71 and tables 2-3). 

According to the results by Merritt et al. (2011), we 
approximate i^emri as: 



^EMRI 




(6) 



where N is the number of stars initially with semi-major 
axis (with respect to the SMBH) a - 2 - 10 x 10~ 3 pc. N 
is basically unknown. We also note that the adopted value 
of iioss was derived for a SMBH-to-BH mass ratio equal 
to 2 x 10 4 . Therefore it is suited for ttismbh = 2 x 10 5 
M Q (adopting m co = 10 M Q ). On the other hand, most of 
the contribution for the detection rates (in Table 1) comes 
from such 'low-mass' SMBHs, as we will discuss in the next 
section. 

We stress that, at the moment of capture, the eccen- 
tricity is likely very high (e.g., Merritt et al. 2011, and ref- 
erences therein; but see Miller et al. 2005 for the effect of 
tidal disruption of binaries), and GWs are emitted in short 
pulses during pericentre passages. Then, the orbit gradually 
shrinks and circularizes (over a timescale of ~ 10 3 — 10 8 yr, 
e.g. Barack & Cutler 2004b), because of GW emission, and 
the GW emission can be observed as continuous. In the next 
Section, we define clso &s the eccentricity at the last sta- 
ble orbit (LSO) and we integrate it back in time accounting 
for GW circularization (see the Appendix A and Barack & 
Cutler 2004a). 

We note also that Merritt et al. (2011) adopt a steady- 
state, cuspy model of the galactic centre (e.g., Sigurdsson 
& Rees 1997; Hopman & Alexander 2006b; Freitag, Amaro- 
Seoane & Kalogera 2006), to derive the estimate that we 
use in equation ([5]). In these models, which are relaxed 
and strongly mass segregated, the mass within ~ 0.1 pc 
of the SMBH is contributed mainly by stellar-mass BHs. 
Therefore, two-body scatterings dominate the orbital evolu- 
tion in proximity of the Schwarzschild barrier (e.g., Freitag 
et al. 2008; Alexander & Hopman 2009), whereas dynamical 
friction is negligible. 

If stars in the central parsec of a galaxy follow a flat 
core rather than a relaxed cusp, the time for a stellar-mass 
BH (initially located out of the core) to reach the cen- 
tre of the host galaxy can easily exceed the Hubble time 
(Merritt 2010). Thus, in flat cores, the number of stel- 
lar BHs within the influence radius of the SMBH is still 
low with respect to the total number of objects, including 
lower-mass remnants and stars. In this case, dynamical fric- 
tion can be more important than two-body scatterings to 
penetrate the Schwarzschild barrier (according to the new, 



2 Background torques are effective only on a timescale shorter 
than the fastest mechanism that changes the relative orientation 
of a star with respect to the gravitational field of the background 
stars. It can be shown that, for high eccentricity, Schwarzschild 
precession is the fastest mechanism (Merritt et al. 2011). 
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more general expression of dynamical friction, derived by 
Antonini & Merritt 2012). This likely leads to lower EMRI 
rates in the case of flat cores and deserves further studies. 

Whether the galaxies hosting a NC have a central stellar 
density cusp or a flat core is an open question. The MW, 
which was long believed to have a steeply-rising mass den- 
sity near the SMBH, was recently found to have a parsec- 
scale core in the old red giant branch (RGB) stellar com- 
ponent (Buchholz, Schodel & Eckart 2009; Do et al. 2009; 
Bartko et al. 2010). Other MW-like galaxies, hosting both 
a NC and a SMBH at their centre, might behave in the 
same way. Therefore, steady-state mass segregated models 
might be inappropriate to describe a MW-like galactic cen- 
tre. On the other hand, only the RGB stars were observed 
to follow a flat core distribution in the central parsec of the 
MW, while the non-relaxed young stars show a rather cuspy 
profile (e.g., Do et al. 2009). We do not know whether the 
old main-sequence stars and especially the stellar remnants 
follow the same cored distribution as RGB stars (e.g., Dale 
et al. 2009; Davies et al. 2011). 

A further possible caveat is that a currently cored pro- 
file might have been cuspy in the past, and then the cusp 
was destroyed as a consequence of dynamical effects (e.g., 
Merritt 2010; Yusef-Zadeh, Bushouse & Wardle 2012). At 
the moment, the most likely scenario to explain the re- 
moval of the cusp and its transformation into a core is that 
of binary scouring as a consequence of a galaxy merger 
(Gualandris & Merritt 2012). In this case, the central cusp 
of stellar BHs reforms in a time that is at least one relax- 
ation time, and the number of BHs in the proximity of the 
SMBH is smaller than in relaxed multi-mass models. 

In summary, whether galaxies hosting NCs have a cen- 
tral stellar cusp or a core, and for how long these profiles 
survive are very uncertain issues. Our predictions for the 
detectable EMRI rate, based on a mass-segregated cuspy 
model, represent the most optimistic case and should be 
regarded as upper limits. 

3. Results 

3.1. Maximum redshift 

Figs. [TJ [5] and [3] show the maximum redshift (z max ) as a 
function of ttismbh in the case of LISA, NGO and DECIGO, 
respectively. In these Figures, to — 5 yr and Xmisskm = 5 yr. 
Here, the value to is the assumed time before the plunge-in, 
i.e. the time elapsed between the beginning of the observa- 
tion with the gravitational interferometer and the moment 
in which the companion reaches the LSO (in this paper, we 
study the cases where to = 5, 6, 10 and 15 yr). T mission is 
defined as the duration of the mission. In this paper, we 
consider Tmission = 5 and 2 yr as two fiducial values of the 
mission duration, for the three considered gravitational ex- 
periments (see, e.g., Gair et al. 2010). We remind that z max 
was obtained for a sky-location and orientation averaged 
signal-to-noise ratio higher than 10 (see Section 2.1 and 
the Appendix A). 

In the case of LISA (Fig. [TJ , z max reaches high values 
(z ~ 2) for msMBH ~ 10 5 — 1O 6 M and eccentricity at 
the LSO eLSO < 0.5. The value of z max for LISA drops 
below ~ 1 for mgMBH < 10 5 M Q and for ttismbh > 2 x 
1O 6 M . The behaviour of z max for NGO (Fig. [2]) is similar 
to that of LISA, in the sense that z max drops quite fast 
for msMBH < 10 5 M Q and for tosmbh > 2 x 10 6 M Q . On the 
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Fig. 1. z max in the case of LISA, as a function of wsmbh for 
three different simulations. In all the shown models, j = 1 
and Tmission = to — 5 yr (see Table 1). Solid line (red on 
the web): model with eLSO = 0.3. Dashed line (blue on the 
web): model with eLso = 0.5. Dotted black line: model with 
eLSO = 0.7. 
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Fig. 2. z max in the case of NGO, as a function of ttismbh for 
three different simulations. In all the shown models, j = 1 
and Tmission = to = 5 yr (see Table 1). Solid line (red on 
the web): model with eLso = 0.3. Dashed line (blue on the 
web): model with eLSO = 0.5. Dotted black line: model with 
eLSO = 0.7. 

other hand, the maximum value of z max for NGO (obtained 
for msMBH ~ 10 5 — 10 6 M Q ) is z max ~ 0.3, much lower than 
in the case of LISA. 

In the case of DECIGO (Fig. [3]), z max ~ 3 for ttismbh ~ 
10 5 M Q , but z max keeps rising for smaller BH masses and 
reaches a value of 8 — 10 for tosmbh ~ 5 x 10 3 M Q . We 
stress that the existence of BHs in the 10 3 — 10 5 M mass 
range is still controversial (e.g., Miller & Colbert 2004, and 
references therein). 

3.2. Detection rates 

Tables 1, 2 and 3 summarize the estimated values of the 
detection rate R for the simulated models, in the case of 
LISA, NGO and DECIGO, respectively. 

Fig. H] shows the behaviour of R as a function of eLso in 
the case of LISA, NGO and DECIGO. In the top panel of 
Fig. |H the ttismbh — TfT S ph relation follows the prescriptions 
by G12a, while in the bottom panel it follows the linear 
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Table 1. Simulated values of the rate R for the old LISA 
configuration. 
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Notes. The tabulated values of R depend on the spin parameter 
(J), on the eccentricity at the last stable orbit (eLSo), on the 
initial time of observations, defined as time before the plunge- 
in (to), on the expected duration of the mission (T m i ss i on ), on 
p = m sp h/mh (see Section 2.1 and Appendix B), and on /bh 
(i.e. the adopted thsmbh —m sp h relation). The value of R out of 
(inside) parentheses refers to /bh from G12a (MH03). Tabulated 
rates are for N = 1 and t\ OBB — 10 7 yr. 



Table 2. The same as Table 1, but for the NGO configu- 
ration. 
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Nominal duration of the NGO mission is 2 yr. 



view of EMRIs in NCs 
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Fig. 3. z max in the case of DECIGO, as a function of 
mgMBH for three different simulations. In all the shown 
models, j = 1 and T m i ss i n = to = 5 yr (see Table 1). Solid 
line (red on the web): model with glso = 0.3. Dashed line 
(blue on the web): model with eLSo = 0.5. Dotted black 
line: model with eLSO = 0.7. 



Table 3. The same as Table 1, but for DECIGO. 
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scaling derived by MH03. Tables 1, 2, 3 and Fig. 0] indicate 
that LISA and DECIGO are expected to detect a factor of 
> 100 more EMRI events than NGO. 

A comparison among Tables 1, 2 and 3 indicates that 
LISA, NGO and DECIGO have similar responses to the 
physical parameters influencing EMRI events (although the 
overall sensitivity of NGO is much lower than that of the 
other two instruments). In particular, for all the considered 
space-borne instruments, R is slightly affected by the spin 
parameter. R mildly depends on the eccentricity (Fig. [4]): 
values of eLSO > 0.7 often correspond to lower values of R, 
since z m ax is lower for er,so > 0.7. 

Tables 1, 2 and 3 show that the detection rate is very 
sensitive to the parameter to, which indicates how far the 
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Fig. 4. The total detection rate R versus the adopted ec- 
centricity at the LSO (eLso), hi the case of LISA (filled cir- 
cles, red on the web), DECIGO (filled triangles, blue on the 
web) and NGO (filled squares). For each case, j = 1, to — 5 

yr, ^mission = 

5 yr (see Tables 1, 2 and 3), p = 1.06 x 10 2 , 
N = 1 and t loss = 10 7 yr. Top panel: m S MBH - m sph rela- 
tion set according to G12a. Bottom panel: tosmbh — TO sp h 
relation set according to MH03. 
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Fig. 5. Cumulative redshift distribution of the expected 
detections per year, in the case of LISA (dotted line, red on 
the web), DECIGO (dashed line, blue on the web) and NGO 
(solid line). For each case, j — 1, glso = 0.1, t — 5 yr, 
T missi o„ = 5 yr (see Tables 1, 2 and 3), p = 1.06 x 10~ 2 , N = 
1 and t\ oss = 10 yr. Top panel: tosmbh — TO sp h relation set 
according to G12a. Bottom panel: tosmbh — TO sp h relation 
set according to MH03. 



observed systems are from final coalescence. R rapidly de- 
creases as to increases. To accurately derive R, we should 
'average' our results over a distribution of values of to - This 
is beyond the aims of this paper. Therefore, we will consider 
the estimates of R obtained for to = 5 yr as upper limits, 
as assuming to — 5 yr maximizes the probability of ob- 
serving a merger event for all the duration of the mission. 



Furthermore, the detection rate sensibly depends on the 
duration of the mission (^mission)- If we assume that LISA, 
NGO and DECIGO will operate only for two years and that 
t — 5 yr, the expected rates are lower by a factor of « 10. 

In our approximate model, R scales linearly with vemri- 
this quantity is highly uncertain, but our estimates can be 
rescaled for different choices of i/emri (the adopted values 
can be considered as upper limits). Finally, the halo-to- 
spheroid mass ratio (p) changes dramatically the rate R, 
and is another poorly constrained factor. 

From Fig. [Hand from Tables 1, 2 and 3, we note also that 
there is a strong dependence of R on the assumed model 
for /BH(TO S ph), that is on the assumed tosmbh — TO sp h re- 
lation. In the case of both LISA and NGO, the relation 
derived by G12a predicts a factor of s=s 10 less detections 
than the linear scaling obtained by MH03. This is not true 
for DECIGO, whose expected detection rate is almost un- 
affected by the assumed /BH("isph)- 

The reason of this behaviour is quite intuitive: the 
tosmbh— TO sp h scaling relation obtained by G12a (tosmbh oc 
TOgph) implies that relatively low mass BHs (< 10 5 M©) are 
associated with most of the spheroids in the 10 s — 10 10 Mq 
mass range, whereas the linear scaling relation derived by 
MH03 predicts that SMBHs with 2 x 10 5 M Q < to S mbh < 
2 x 10 7 M© are hosted by spheroids in the same mass range 
as above. LISA and NGO have very low sensitivity to 
EMRIs with tosmbh < 10 5 M© (Figs. [TJ and HJ), and this 
explains the low values of R when the relation by G12a 
is assumed. Instead, DECIGO is particularly suitable to 
observe EMRIs involving 10 4 M BHs (Fig.©, while per- 
forming as well as LISA for tosmbh ~ 10 5 - 10 6 M Q . This 
explains why the predicted detection rate for DECIGO is 
almost unaffected by the adopted tosmbh — TO sp h scaling 
relation. 

Fig. [S] shows the cumulative redshift distribution of the 
expected detections per year, that is it shows how many de- 
tections per year are expected for EMRIs at redshift < z, 
as a function of z. Fig. [5] confirms that NGO is blind for 
EMRIs with z > 0.3 (as a consequence of the behaviour of 
Zmaxi see Fig. From Fig. [5] we note also the very differ- 
ent behaviour of LISA and DECIGO with respect to the 
adopted tosmbh — TO sp h scaling relation. When the MH03 
relation is assumed, LISA and DECIGO have almost the 
same behaviour, in terms of both total detection rate and 
dependence on redshift: most of detectable EMRIs occur at 
z ^ 2. Instead, when the G12a relation is assumed, most of 
detectable EMRIs occur at z £ 2 and z £ 4 for LISA and 
DECIGO, respectively. We stress that our predictions based 
on the Press-Schechter formalism might overlook many im- 
portant processes affecting BH and galaxy evolution at 
z » 1. 



3.3. Comparison with previous studies 

In this section, we compare the rates derived for the orig- 
inal LISA configuration and for NGO with those found 
in previous studies (e.g., Gair et al. 2004; Rubbo, Holley- 
Bockelmann & Finn 2006; Gair et al. 2010; Amaro-Seoane 
et al. 2012a). In particular, we will focus on Gair et al. 
(2004) and on Amaro-Seoane et al. (2012a), for LISA and 
NGO, respectively. 

Gair et al. (2004) find a EMRI detection rate - 30-350 
yr~ x (for a 10 M© BH), depending on the assumed prop- 
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erties of the LISA mission. This result is quite consistent 
with our predictions, when we assume the G12a model, and 
is a factor of ~ 10 lower than our results, when we adopt 
the MH03 prescriptions. On the other hand, there are sig- 
nificant differences between our method and that by Gair 
et al. (2004). 

Firstly, Gair et al. (2004) assume a EMRI rate i'emri ~ 
10~ 6 (m SM BH/3 x 10 6 ) 3 / 8 yr" 1 per SMBH (for mergers with 
a 10 Mq BH). This is a factor of 10 higher than our assumed 
EMRI rate (see equation [5]) , when we adopt N — 1 (as in 
our Figs. H] and [5]). 

Secondly, Gair et al. (2004) assume a different mass 
function for the SMBHs with respect to us. In particular, 
Gair et al. (2004) estimate j^mbh from the msMBH _ a 

relation (where a is the velocity dispersion, e.g., Ferrarese 
& Merritt 2000; Gebhardt et al. 2000; Graham et al. 
2011) and from the luminosity versus u relation (Aller & 
Richstone 2002). Instead, we use the msMBH — "^sph rela- 
tion (e.g., G12a) and the Press-Schechter formalism (Press 
& Schechter 1974). Furthermore, we consider only those 
SMBHs that can coexist with a NC (otherwise the cap- 
ture rate would likely be much lower), whereas Gair et al. 
(2004) take into account all the BHs in the mass range 
10 5 — 5 x 10 6 M Q , excluding only those that are located in 
Sc and Sd galaxies. 

These differences are important, in terms of EMRI rates, 
as it can be shown by calculating the expected merger rate 
per unit volume at redshift z = (^emri = ismbh(z = 
0) i-'emrI) where iismbh(^ = 0) is the number density of 
SMBHs at z = for the considered mass range and i^emri 
is the EMRI rate per SMBH). In our model, considering 
only SMBHs with 10 5 < m SM BH < 5 x 10 6 M (i.e., the 
mass range adopted by Gair et al. 2004), j>emri ~ 20 Gpc~ 3 
yr _1 and #emri ~ 2 Gpc~ 3 yr^ 1 , assuming the MH03 and 
the G12a scaling relation, respectively. 

Our estimate of j^emri derived from the G12a scaling 
relation is similar to the one reported in table 2 of Gair et 
al. (2004), while it is one order of magnitude higher than 
the one reported by Gair et al. (2004) when the MH03 
scaling relation is adopted instead. Our estimate of Cemri 
based on G12a, although very similar to the one derived 
by Gair et al. (2004), comes from the product between our 
merger rate ^emri, which is a factor of 10 lower than the 
one adopted by Gair et al. (2004), and the number density 
of SMBHs with 10 5 < m SM BH < 5 x 10 6 M Q , which is 
a factor of 10 higher than the one in Gair et al. (2004). 
We stress that our value of £emri is an upper limit, as 
^emri was obtained assuming a cuspy density profile inside 
the SMBH influence radius (see the discussion in Section 
2.2) and we impose that /bb = 1, i.e. that all NCs in the 
considered mass range host a SMBH. 

In addition, Gair et al. (2004) do not make a fully 
consistent cosmological integration in a ACDM Universe. 
Furthermore, the number of detectable events reported in 
table 3 of Gair et al. (2004) was derived imposing that 
z max cannot be higher than one. On the other hand, Gair 
et al. (2004) use the same kludge waveforms we adopt (from 
Barack & Cutler 2004a) and include an accurate treat- 
ment for detection statistics (while we consider an averaged 
signal-to- noise ratio over the duration of the mission). To 
check the importance of the redshift integration, we can 



substitute our equation ((5]) with a less accurate calculation 
(see the appendix C of Mapelli et al. 2010): 

R = "smbh(z = 0) i^EMRI V c (z cut ), (7) 

where ?^smbh {z = 0) is the density of SMBHs surrounded 
by NCs at redshift z = in the 10 5 - 6 x 10 6 M Q range, 
^emri is defined in Section 2.2, and V c (z cut ) is the comov- 
ing volume up to redshift z cut . For comparison with Gair et 
al. (2004), we assume that z cu t is equal to the minimum be- 
tween 1 and Zmax- From equation (J7J, we obtain an approxi- 
mate EMRI detection rate R « 2 x 10 3 yr" 1 and R w 2 x 10 2 
yr _1 , for the MH03 and the G12a scaling relation, respec- 
tively. These values are of the same order of magnitude as 
those listed in our Table 1, indicating that the accurate 
redshift integration does not affect significantly the rate in 
the case of LISA. 

We now consider the predictions for NGO, by comparing 
our results with those from Amaro-Seoane et al. (2012a). 
In their fig. 21, Amaro-Seoane et al. (2012a) show z max as 
a function of the SMBH mass, for EMRIs with a 10 M 
BH. Our Fig. [2] is in fair agreement with the sky-averaged 
horizons shown in fig. 21 of Amaro-Seoane et al. (2012a), 
taking into account that Amaro-Seoane et al. (2012a) sky- 
averaged horizons were calculated for Teukolsky waveform 
models (Teukolsky 1973), T m i ssion = 2 yr, random gener- 
ated plunge times < £o/ vr < 5 and (SNR) > 20, whereas 
in our Fig.[2]we adopt Barack & Cutler (2004a) waveforms, 
T m i SS ion = = 5 yr and (SNR) > 10. The expected EMRI 
detection rates for NGO, that can be derived from table 2 
of Amaro-Seoane et al. (2012a) (R « 15 — 25 yr" 1 ) are con- 
sistent with the values reported in our Table 2, but quite 
optimistic for a 2-yr mission. 

4. Conclusions 

NCs seem to be common in galaxies, independently of the 
Hubble type. Recent studies (Seth et al. 2008a, 2008b; 
Gonzalez Delgado et al. 2008, 2009; Graham & Spitler 2009) 
indicate that NCs co-exist with SMBHs in the same nu- 
cleus. The fraction of galaxies hosting both a NC and a 
SMBH is particularly high among spheroids with stellar 
mass ranging from ~ 10 8 M Q up to a few ~ 10 10 M Q . Given 
their high stellar densities and high escape velocities, NCs 
are expected to retain a large fraction of their stellar-mass 
BHs. Thus, it is reasonable to expect that the rate of EMRIs 
is enhanced in NCs. 

We consider sensitivity curves for three different 
space-based GW laser interferometric mission concepts: 
the Laser Interferometer Space Antenna (LISA), the 
New Gravitational wave Observatory (NGO) and the 
DECi-hertz Interferometer Gravitational wave Observatory 
(DECIGO). 

We show that GW signal from SMBH— BH mergers can 
be observed as far as redshift z ~ 2 and z ~ 8, for LISA 
and DECIGO, respectively, while this limit is substantially 
lower (z ~ 0.3) for NGO. We predict that, under the most 
optimistic assumptions, LISA and DECIGO will detect up 
to thousands of EMRIs in NCs per year, while NGO will 
observe up to tens of EMRIs per year. 

The estimate of the detection rate R for SMBH— BH 
EMRIs is subject to a plethora of severe uncertainties. 
First, the EMRI rate depends on the dynamics and on the 
relativistic effects in the neighborhoods ( 10~ 2 pc) of the 
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SMBH. In particular, the collisional dynamics of relativistic 
star clusters around SMBHs is poorly understood. Our pre- 
dicted rates are upper limits, in the sense that the EMRI 
rate was derived assuming a mass-segregated cusp in the 
proximity of the SMBH (Merritt et al. 2011). The exis- 
tence of a flat core rather than a mass-segregated cusp in 
the central parsec of galaxies can severely affect the EMRI 
rate (e.g., Antonini & Merritt 2012). 

Furthermore, the cosmological rate of EMRIs depends 
on the dark-to-luminous mass ratio in galaxies, which is 
highly uncertain. In addition, we assumed that the fraction 
of spheroids that host both a SMBH and a NC is fbb = 1 
for a spheroid mass range 10 s < m sp h/M0 < 10 10 . This 
is consistent with recent observations (Graham & Spitlcr 
2009), but is still based on a small sample of galaxies (~ 
20 — 30 galaxies are known to host both a SMBH and a 
NC in their nucleus). If the available sample suffers from 
any biases, the actual fbb may be < 1, although values of 
fbb << 1/2 are unlikely. In addition, /bb might change with 
redshift. In general, we do not know how NCs and SMBHs 
evolve with redshift and we can only assume that the halo 
occupation of NCs and SMBHs is the same up to z ~ 2. 

The mass distribution of SMBHs is poorly understood, 
especially in low- mass galaxies. Previous studies (e.g., Gair 
et al. 2004; Amaro-Seoane et al. 2012a) derive the SMBH 
mass distribution from the tosmbh — c relation, whereas 
we use the tosmbh — m-sph relation. Both relations are still 
largely uncertain at the low-mass tail. In the case of the 
tosmbh — "i sp ii relation, while previous studies proposed 
a linear scaling (e.g., MH03; Haring & Rix 2004), more re- 
cent data indicate a steeper scaling (ttjsmbh oc to^, G12a, 
G12b). This difference implies large discrepancies (by a fac- 
tor of 10) in the expected GW detection rate. 

Finally, the sensitivity curves of space-borne interfer- 
ometers will likely be modified in the next few years, as 
DECIGO is still in the design phase, while the LISA mis- 
sion concept was recently abandoned and the NGO mission 
concept was not selected by the European Space Agency 
(ESA). Differences in the sensitivity curves might strongly 
affect the predictions that we have computed here. 

Therefore, to give an accurate estimate of the detection 
rate R of EMRIs in NCs is beyond the scope of this paper. 
On the other hand, our results in Tables 1, 2 and 3 and 
in Figs. S] and [5] can be considered as upper limits, and 
indicate that EMRIs in NCs might be detected by space- 
born detectors. In summary, the co-existence of SMBHs and 
NCs can boost significantly our chance of observing GWs 
from EMRIs in the galactic centres, although more studies 
of the dynamics of relativistic star clusters are needed, to 
quantify this with accuracy. 
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Appendix A: Method to estimate z max 

We define z max (wsmbh ? m co ) as the maximum redshift at which an event can be detected with a sky location and 
orientation averaged signal-to-noise ratio (SNR) > 10 by a single interferometer. Mergers of SMBHs with stellar BHs 
are classified as extreme mass-ratio inspirals (EMRIs). For such systems, one can estimate (SNR) following the approach 
and waveform approximation described by Barack & Cutler (2004a). The angle-averaged signal-to- noise ratio (SNR) for 
a single synthetic two-arm Michclson detector is given by: 

<^>=£/^f^/,, (t >>, (A,, 



where f n (t) indicates the frequency of the different harmonics at multiple integers of the orbital frequency (in our 
calculations, we adopt n — 1, 20, which safely accounts for all the GW power) and is defined to be f n (t) = nv(t)+j(t)/TT 
(where v{t) is the orbital frequency and j(t) is the direction of pericentre with respect to L x S -i.e. the vector product of 
the orbital angular momentum unit vector L and the SMBH spin direction S). h c , n (f n (t)) is the characteristic amplitude 
and has been derived following equation (56) of Barack & Cutler (2004a; see also Cutler et al. 1993, 1994). Finally, 
Sh{fn{t)) accounts for the total noise of the interferometer and we give expressions for LISA, NGO and DECIGO in the 
next subsection. In this Paper, we derive f n (t) by solving, backwards in time, the system of differential post- Newtonian 
equations (27-30) of Barack & Cutler (2004a). In particular, we evolve the mean anomaly $(i), the orbital frequency 
v{t), the direction of pericentre j(t) and the orbital eccentricity e(t) starting from the LSO, backwards in time down 
to the assumed starting time of the observation. Wc define as io the duration elapsed between the starting time of the 
observation and the epoch when the companion reaches the LSO (in this paper, we consider four cases: t = 5, 6, 10 
and 15 yr). We furthermore define T m ; ss i on the duration of the mission (in this paper, we consider two cases: T m j ss i on = 5 
and 2 yr). As arbitrary initial conditions, we set $lso = 0, 7lso = 0> and we keep the LSO eccentricity eLso as a free 

parameter (see Table 1), which yields ^lso = (2 ttwsmbh) 1 [(l — e Lso) / (6 + 2 eLso)] 3 ^ 2 - 

Using the above equations, we can derive the luminosity distance Dl (z max (msMBH, m co )) at which an event can be 
detected with a SNR > 10 by a certain interferometer (see, e.g., Mapelli et al. 2010 and references therein). We can, 
thus, derive z max (msMBH, m co ) by inverting the expression of the luminosity distance in a fiat ACDM model: 



-DL(z max (TOSMBH, m c°)) 



Tr(l + z max(mSMBH, ™co)) / T^T • (A.2) 

■"0 Jo t\ z ) 



A.l. LISA noise 

Sh{fn(tj) accounts for the total noise of LISA noise and is defined as (equation 54 of Barack & Cutler 2004a): 
S h (f n (t)) = sT t+ ^(f n (t)) + Sr sa '(/n(i)), (A.3) 

where S]^ st+sal (f n (t)) is defined in equation (52) of Barack & Cutler (2004a, see also Hughes 2002) and accounts for 

both the instrumental noise (Sj^ st (f n (t)) , Finn & Thornc 2000) and the Galactic confusion noise (Sf^ (f n (t)), Nelemans, 
Yungelson & Portegies Zwart 2001, see also Bender & Hils 1997): 



sr t+gal (/«w) 



5rt(/ " (t)) — sT\.ut)) + sf(f n (t)) 



i . ,, h \ j „ v , , ■ 7 i j „ w . (A. 4) 

_exp(- ^T mission N f ) 

where n ~ 4.5 (Cornish 2003), T mission yr is the duration of the mission, Nf = 2 x 10 _3 Hz _1 (Hz/.f„(t)) n / 3 (Hughes 

2002). S^ st (f n (t)) is the same as in equation (48) of Barack & Cutler (2004a). 5^ al (/„(i)) is the same as in equation 

(51) of Barack & Cutler (2004a). Finally, S^ gal (f n (t)) accounts for the extragalactic white dwarf background (Farmer 
& Phinney 2003) and is the same as in equation (50) of Barack & Cutler (2004a). 



A.2. NGO noise 

In the case of NGO, S™ st (f n (t)) is defined as (equation 5 of Amaro-Seoane et al. 2012a): 

fn(t) 



Cinst^ / i \ \ _ 20 4S' x , acc (/„(£)) + S'x,sn(/n(i)) + Sx^omn (fn (t)) , 
&h [Jn{t)) - "Tj J2 \ 1 + 



0-41 (2%) 



(A.5) 



where L = 10 9 m, S^ acc {f n {t)) = 1.37 x 10~ 32 (l + {j^f m 2 Hz" 1 is the power spectral density of the 

residual acceleration of the test mass, S XtSn (f n (t)) = 5.25 x 10~ 23 m 2 Hz -1 is the shot noise power spectral density, and 
Sx,omn{.fn{t j) = 6.28 x 10~ 23 m 2 Hz -1 is the power spectral density of the other possible measurement noises. 

As in the case of LISA, we combine the instrumental noise S™ at (f n (t)) with the Galactic confusion noise and with 
the extragalactic white dwarf background (see Section A.l). On the other hand, the Galactic confusion noise and the 
extragalactic white dwarf background are negligible compared to S™ st (f n (t)) in the case of NGO. 
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A.3. DECIGO noise 

In the case of DECIGO, we calculate the noise as in equation (1) of Yagi & Tanaka (2010): 



Sh(fn(t)) = min 



Sr\fu(t)) 



S e * sal (f n (t))K(f n (t))+0M Sr(Ut)), (A.6) 



.exp(-KT mi 1 ssion 7V / ) 

where the instrumental noise spectral density (S l £ at (f n (t))) is given by (Yagi & Tanaka 2010) 



, .,r,,u .r., ,,, is fl , r 2, . 2-3 x 10-7 2.6 x 10- 81 

1 ; x 4 (l+x 2 ) a: 4 



ST(f n (t)) = 5.3 x 10- 



Hz" 1 , (A.7) 



where x = f n (t)/f p with f p = 7.36 Hz. k, Tmissiom Nf, Sf^ (f n (t)) and 5 , ° xga (/„(£)) are the same as defined for LISA 
in Section fA. 11 Finally, the factor TZ(f n (t)) = exp [—2 (/„(t)/0.05 Hz) 2 ] represents the cutoff of the white dwarf— white 
dwarf binary confusion noise (Yagi & Tanaka 2010). 

Appendix B: Estimate of p from the Press-Schechter formalism 

We can derive p in the following way. Driver et al. (2007) estimate that the current stellar mass density of spheroids is 
Psph ~ 3.0 x 10 s hM Q Mpc -3 (where h = 0.71, Larson et al. 2011). Using the Schechter function (see Driver et al. 2007 
for details), we find that the mass fraction of spheroids with mass 10 8 < m sp h/M0 < 10 10 is g ~ 0.34. Thus, the current 
mass density of halos hosting spheroids with mass 10 8 < m sp h/M0 < 10 10 is 



f* \ J \ Psph + Pdisc 

where /» is the fraction of baryons in stars, Pdisc is the current stellar mass density of discs, J7m = 0.27 and flh — 0.044 
(Larson et al. 2011). Assuming g = 0.34, /* = 0.119 h and p disc = 4.4 x 10 8 h M Mpc -3 (from Driver et al. 2007), we 
find p h = 1.3 x 10 10 M Mpc" 3 . 

A density — 1.3 x 10 10 M Q Mpc -3 can be obtained by integrating the Press-Schechter function (Press & Schechter 
1974; Eisenstein & Hu 1998, 1999) between m h i = 9.44 x 10 9 M Q and m h2 = 9.44 x 10 11 M Q (which satisfies the further 
requirement that = lOOrrihi)- Whi = 9.44 x 10 9 M Q and = 9.44 x lO n M0 can be used as integration limits in 
equation ([5]). Finally, we can derive also the value of p as p = m sp hi/mhi = 1.06 x 10~ 2 . 

We note that this value of p matches very well with the observations of the Milky Way, but there may be large 
deviations from it, depending on the morphological type (see, e.g., Graham & Worley 2008). A different value of p may 
significantly affect the results (see the examples in Tables 1, 2 and 3). 
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